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Guderley  acting  as  project  engineer. 
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ABSTRACT 


The  wall  pressure  distribution  in  a choked  tunnel  is  computed, 
following  a method  outlined  previously.  By  the  calculus  of  residues 
the  solution  is  brought  into  such  a form  that  its  structure  is  displayed 
more  clearly. 
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INTRODUCTION 


For  experiments  on  the  wall  influence  in  transonic  flow  the  wall 
pressure  distribution  in  a choked  tunnel  has  sane  interest.  It  is 
determined  in  the  present  report  under  the  assumption  that  the  wind 
tunnel  is  very  wide  in  comparison  to  the  model  dimensions. 

The  paper  is  closely  related  to  an  investigation  by  B.  Marschner 
(Reference  1)  on  the  influence  of  the  tunnel  wall  on  the  pressure 
distribution  over  the  body.  Indeed  his  analysis  contains  in  principle 
all  ideas  necessary  for  the  determination  of  the  wall  pressure  distri- 
bution. However,  the  solution,  thus  obtained,  would  appear  in  a form 
which  is  not  yet  suited  for  the  numerical  evaluation  and  does  not  put 
into  evidence  its  essential  properties.  The  main  task  of  the  present 
report  is  a transformation  of  the  original  result  into  a more  suitable 
form  by  means  of  the  calculus  of  residues. 

An  aerodynamicist  will  consider  as  essential  in  this  paper  the 
pressure  distribution  over  the  wind-tunnel  wall.  But  the  mathematical 
part  is  also  quite  enjoyable  because  of  the  spectacular  simplifications 
which  are  achieved. 

Since  the  problem  is  explained  in  full  detail  in  the  paper  of 
Marschner,  and  since  it  follows  the  procedure  outlined  in  Reference  2, 
the  general  explanations  can  be  kept  rather  brief. 
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SECTION  I 


THE  BOUNDARY  VALUE  PROBLEM 


Figure  I shows  the  flow  in  the  physical  plane.  If  the  tunnel  is 
choked,  a sonic  line  (OC)  will  extend  frcm  the  body  to  the  wind-tunnel 
wall.  Mach  waves  which  start  from  the  supersonic  contour  of  the  body 
will  either  end  at  the  sonic  line  or  at  the  wind-tunnel  wall.  The  one 
Mach  wave  which  separates  these  two  classes  is  called  the  limiting 
Mach  wave,  since  it  gives  the  downstream  limit  of  the  part  of  the 
supersonic  flow  field  which  influences  the  subsonic  flow  (line  OD). 

The  hodograph  picture  is  shown  in  Figure  2.  To  fix  the  ideas  the 
body  has  been  assumed  to  be  a wedge.  Corresponding  points  in  the 
physical  plane  and  in  the  hodograph  plane  have  the  same  notation.  The 
surface  of  the  wedge  maps  into  the  line  BC,  the  wind-tunnel  wall  into 
the  line  AO  (with  different  values  of  the  stream  function  on  the  upper 
and  lower  side  of  this  line,  since  they  correspond  to  the  upper  and 
lower  walls  of  the  tunnel ) • At  point  A a singularity  arises  since  it 
corresponds  to  the  flow  in  the  tunnel  at  >a  large  distance  upstream  and 
thus  all  streamlines  originate  at  this  point.  AB  corresponds  to  the 
line  of  symmetry  of  the  wind  tunnel;  along  this  line  the  value  of  the 
stream  function  is  zero. 

The  sharp  corner  which  occurs  at  C in  the  physical  plane  maps  into 
the  characteristic  CD,  which  will  be  referred  to  as  the  shoulder 
characteristic. 

In  the  transonic  approximation  one  usually  introduces  instead  of 
the  coordinate  w (w  absolute  value  of  the  velocity  vector)  another 
variable  £ , defined  by 

_ 1 

w » w*  (1  + ( t + 1)  ^ £ ) (l) 

where  w*  is  the  sonic  velocity  and  y is  the  ratio  of  the  specific 
heats.  The  other  independent  variable  is  S , the  angle  of  the  velocity 
vector  with  the  horizontal  axis.  Actually,  a limiting  process  is 
carried  out  in  which  it  is  assumed  that  the  body  becomes  very  slender 
and  the  Mach  number  very  close  to  one.  By  this  limiting  process  the 
value  of  £ for  the  stagnation  point  moves  to  - co.  The  boundary  value 
problem  in  the  £ , & -plane,  is  shown  in  Figure  3« 

This  problem  will  be  treated  with  a method  described  in  Reference 
2.  The  approach  chosen  has  the  advantage  that  the  solutions  appear  in 
a particularly  simple  form  in  the  limiting  case,  where  the  choking  Mach 
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number  is  close  to  one;  or,  in  other  words,  if  the  wind-tunnel  dimen- 
sions are  large  in  comparison  to  the  model* 


SECTION  II 

PARTICULAR  SOLUTIONS 


In  the  transonic  approximation  the  equation  for  the  stream 
function  y assumes  the  form 


= o 


(2) 


Introducing  new  coordinates 


t 


(iej 


Vj 


- 1*  + a g) 


(3a) 


(3b) 


one  obtains  as  the  equation  for  the  stream  function 


■ fff'-p/f  +iy  ’ fr?rr  * *&y 


(4) 


Particular  solutions  of  the  form 


y * G(}) 


(5) 


have  been  studied  in  Reference  2*  They  yield  the  differential 
equation  for  G 

- i £(>-?) S'  + (-*  *i)f  = o 

/ (6) 
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■which  can  be  Witten  as 


fyo-tr-etj 


0-?) 


VC 


G=o 


(7) 


The  function  g(  ^ ) is  expressed  by 

■/i  - j ^ 

$(<!>  ’ f 

The  function  G depends  upon  / and  j therefore,  the  notation 
?<  L * l mil  be  used  occasionally.  The  arising  particular  solutions 
for  y>  -will  be  either  symmetric  or  antisymmetric  with  respect  to  the 
axis  0 ■ 0;  therefore  superscripts  s,  symmetric,  and  a,  antisymmetric, 
will  be  used. 


Let  C2  be  a constant  whose  choice  will  be  discussed  later.  If 

along  the  lines  £ * - co  and  J **  Cg  the  boundary  condition  y>  = 0, 

and  consequently  G = 0,  is  prescribed,  the  f -unctions  G and  the  para- 
meters A are  determined  by  an  eigen  value  problem.  It  deviates  from 
classical  eigen  value  problems  only  by  the  fact  that  the  coefficient  of 
G in  Equation  (6)  changes  its  sign,  namely,  at  f » 0,  i.e.,  with  the 
transition  from  the  subsonic  to  the  supersonic  range.  One  obtains  an 
infinite  number  of  positive  and  an  infinite  number  of  negative  eigen 
values  and  a corresponding  system  of  eigen  functions  G.  The  eigen  values 
may  be  arranged  according  to  their  magnitude  and  denoted  by 


^ ....  1 ^ ....  A .... 

“ - j -JL  - / / -2  J H 

where  the  negative  subscripts  refer  to  negative  eigen  values  and  the 
positive  subscripts  to  positive  eigen  values.  The  eigen  functions  G 
may  be  denoted  accordingly. 

The  system  of  eigen  functions  can  be  shown  to  be  complete.  One 
has  the  following  relations  of  orthogonality 


0-}3)w 


Gh  ■ Gk 


O 


for  h / k 


(2) 


The  particular  solutions  for  y>  which  can  be  formed  with  these 
expressions  G will  be  denoted  in  the  following  manner: 
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(9a) 


Yu  * / • 6h 


hh  m ? • Gh 


It*  - f * ^ (3  fcjjyi) 

Yvt,  - 


(9b) 


(9c) 

(9d) 


Here  h is  a positive  integral  number,  ^ is  a suitable  constant. 


Because  of  the  completeness  of  the  system  of  eigen 
function  y/  which  fulfills  for  f < f < ( T 1 311(1  f 


functions,  any 
2 arbitrary- 


constants)  the  differential  equation  (h)  and  the  boimdary  conditions 
y-  » 0 along  J m - co  and  f *=  Cg,  can  be  represented  in  this  region 


by  a superposition  of  expressions  (9).  (See  Reference  2,  Appendix  1.) 


SECTION  III 

THE  CHOICE  OF  THE  VALUE  OF  C2 

Since  in  the  original  boundary  value  problem  no  condition  is 
prescribed  along  the  limiting  characteristic  J ■ 1,  the  application 

of  these  particular  solutions  to  the  present  problem  is  not  quite 
self-evident. 

We  proceed  as  follows:  along  the  shoulder  characteristic  the 

boundary  condition  JU  « 0 is  prescribed.  We  now  define  that  y,  s 0 
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in  the  region  CDG  of  Figure  3*  Then  the  differential  equation  for 
is  fulfilled  everywhere  in  the  region  OABCGDO,  even  at  the 
shoulder  characteristic  CD,  (That  along  this  line  a jump  of  the  first 
derivative  of  occurs,  is  not  in  contradiction  to  the  fulfillment 
of  the  differential  equation,) 

In  order  to  prepare  the  use  of  the  eigen  functions,  we  draw  a 
line  £ “ Cg  (0  ^ C2  < 1)  and  prescribe  (in  contradiction  with  the 

requirements  of  the  boundary  value  problem),  that  «*  0 along 
£ * Cg»  It  has  been  made  plausible  in  Reference  2 that  in  the  limit- 
ing case,  C2  lj the  solution  of  the  boundary  value  problem  with  the 

line  { = C2  as  additional  boundary,  tends  to  the  solution  of  the 

original  problem.  For  a purely  supersonic  problem  this  state  of 
affairs  is  quite  obvious  (see  Reference  2);  the  reasoning  for  the 
mixed  subsonic  supersonic  case  is  less  direct.  For  the  present  problem 
this  is  unessential,  since  the  final  form  of  the  solution  will  justify 
the  procedure.  Actually,  the  present  example  may  serve  as  an  illus- 
tration for  the  limiting  process  02"*’  1* 

Perhaps  same  remarks  regarding  the  physical  nature  of  this 
limiting  process  will  be  welcome.  It  seems  as  if,  by  the  limiting 
process  that  Cg-^1,  the  boundary  condition  y>  ■ 0 is  imposed  not  only 

along  the  shoulder  characteristic  CD  but  also  along  the  limiting 
characteristic  DO,  This  is  in  contradiction  with  the  requirements  of 
Triccmi’s  existence  theorem  where  no  condition  can  be  imposed  along  one 
of  the  characteristics  which  form  the  boundary  of  the  supersonic  region. 
On  the  other  hand  the  boundary  condition  has  a physical  significance. 
Let  us  assume  that  the  flow  field  in  a closed  tunnel  has  been  found  up 
to  the  limiting  Mach  wave.  Now  we  construct  the  rear  part  of  the  body 
in  such  a manner  that  it  does  not  reflect  any  Mach  waves.  Then  the 
region  DQOH  of  Figure  1+  maps  into  the  one  characteristic  DO  of  the 
hodograph.  Thus  the  -surface  in  the  hodograph  plane  will  drop  along 
the  line  DO  suddenly  from  the  values  which  correspond  to  the  limiting 
characteristic,  to  the  value  of  ^ » 0 which  corresponds  to  the  line 
DH  in  Figure  lu  This  sudden  drop  in  the  hodograph  is  the  mapping  of 
the  so-called  "lost"  solutions  ("lost"  if  only  such  hodograph  expres- 
sions are  considered  which  are  continuous  in  w and  & ),  Thus,  the 
requirement  p ■ 0 along  the  limiting  characteristic  has  a physical 
interpretation,  although  it  is  meaningless  for  the  mathematical  formu- 
lation of  the  problem. 
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Also  the  boundary  condition  ^ ■ 0 along  f » Cg  can  be  given 

a physical  meaning,  it  requires  that  the  Meyer  expansion  at  the  shoulder 
is  carried  out  only  up  to  point  K and  that  from  there  on,  a certain 
relation  between  the  pressure  and  the  streamline  direction  is  prescribed 
along  the  zero  streamline.  If  C2  -*  1 the  Meyer  expansion  is  carried 

out  further  and  further,  and  the  condition  along  the  zero  streamline 
comes  closer  and  closer  to  the  condition  assumed  for  the  rear  part  of 
the  body  in  Figure  Thus  the  limiting  process  c2  -*•  1 appears  to  be 
justified  by  a physical  consideration. 


SECTION  IV 

REPRESENTATION  OF  AN  ARBITRARY  FUNCTION  BY  A 
SUPERPOSITION  OF  EIGEN  FUNCTIONS  G 


An  arbitrary  function  f(^)  can  be  represented  by  a Fourier 

series  in  the  eigen  functions  G.  Using  the  relations  of  orthogonality 
(S),  one  finds 


hmi  -°° 

~ A 


hm/  \t,  -00 

where  the  constants  C^  are  given  by 


£ 9 J 


(10) 


Oi) 


In  Equation  (10)  the  limiting  process  that  C2  1 will  be  carried 

out*  Valid  for  any  value  of  A , we  have  the  formula  (Appendix  h> 
Reference  2) 
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T 


dt  m 


m My 

df  di/ 


Here  the  functions  G are  defined  by  the  differential  equation  (7)  and 
their  development  at  the  point  £ ® - oo  which  starts  with  the  term 

(-  | The  functions  G are  expressed  by  hyper  geometric  series  in 

Equations  (IS)  of  Reference  2.  For  the  evaluation  of  the  one  needs 
a representation  of  the  in  the  vicinity  of  £ * 1#  According  to 

Equation  (ISd),  Reference  2,  one  has 


(13) 

If  i-i  the  hypergeometric  series  can  be  replaced  by  1#  The  determi- 
nation of  the  positive  eigen  values  has  been  shown  in  Reference  2*  The 
main  idea  is  that  the  first  term  in  Equation  (13)  tends,  for  real 
positive  values  of  y^T  , to  infinity  if  f -*■  1.  Thus  the  boundary 

condition  G ■ 0 requires  that  the  coefficient  of  the  first  term 
vanish.  Hence 

zf^+i)  - o 
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*T\ 


and 


x.  +M) 


A-  /,  X • 


(lit) 


Thus  the  functions  will  be  given  for  positive  values  of  h as 


G-G  7f,  4 'ftf) 


For  these  values  of  X the  limiting  process  that  C2  -►  1 can  be 

carried  out  in  Equation  (12)  without  difficulty*  Correcting  an  error 
in  Reference  2 (Equation  (A26)),  one  finds 

/jl  7// . ±iz^±j  dz  a 


or 


For  the  negative  eigen  values  one  proceeds  as  follows: 
Introducing 


r /T  . 


LP 


or 


and 


X 

”-4 


f/-  - * 


6 


(15) 


(16) 


(17) 
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■we  -write  Equation  (13)  for  G in  the  vicinity  of  f **  1 as 


6. W I**’ 

u 1 l rfcv  +£) . -fj) 


r(-L»+£)  .rt-cv+i) 


71 


+ c»  i 


For  brevity  we  set 


<t(v} 


rfi)  . A()>)  . e ‘ 

rfi»i/.r a»*i) 

where  A(  y ) and  0~  ( V ) are  real.  Then, 


6!  - z Am  f'-f  ’)*  • ^ A-/V 

At  / “ c2  bbe  boundary  condition  G « 0 yields 

«*/*/*;  - $ *7*)  - ° 


(13) 


(19) 


(20) 


Using  Equations  (17)  and  (20)  and  carrying  out  the  limiting  process 
Cg  -*  1 as  far  as  practical,  one  obtains  from  Equation  (19) 

* s - j '**•*» -p*  Jbc) ’fi' -'ey£J 


- n 


dG  w - 6a  £ u ) 


Because  of  Equation  (20)  the  absolute  value  of  the  sine  function 
is  1,  Thus  the  values  of  C . are  given  by 
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(21) 


C-h  tends  to  positive  infinity  as  £ ->  0.  The  values  of  are 

determined  from  Equation  (20)*  One  finds 


°'(vh) 


- » *■(£ 


For  small  values  of  <f  (02  close  to  1),  the  quantity  will  change 

by  only  a small  amount  if  h changes  by  unity.  Denoting  the  change  of 
V between  two  consecutive  values  of  h as  A V,  , one  thus  finds  from 
the  last  equation 

AVh(<r'-  -+  f as  £ * o 


With  this  formula  the  second  expression  occurring  in  Equation  (10)  can 
be  changed  from  a summation  over  h to  a summation  over  the 

DO 


Z pv.j/a/- L.n  ■ 


T AV*  ' AV  ///r, L . £Y t -$»x)ctr 


Hence  one  finds 
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Thus  Equation  (10 ) yields  finally 


Z J^p7(  ■ G(t  H+Mjyar 


f't) 


IT  I (//*>!&* 


Jt  / A*(»> 

y*  o 


SECTION  V 

THE  SOLUTION  OF  THE  BOUNDARY  VALUE  PROBLEM 


By  a line  g • constant  *»  f Q through  point  A of  Figure  3>  the 

region  in  consideration  is  divided  into  an  inner  portion  ( $ < <?0  ) and 

an  outer  portion  ( g > ?0  )•  For  these  regions  different  representa- 
tions of  the  solution  containing  undetermined  coefficients  will  be 
chosen*  The  coefficients  will  then  be  found  by  matching  the  solutions 
along  the  line  g ■ Q * 

In  the  inner  region  one  has  the  boundary  condition  y « 1 along 
OA,  i.e*,  for  / ■ - oo*  A solution  of  the  differential  equation  (1+) 
which  fulfills  this  boundary  condition  is  obtained  for  A » l/l6,  Then 
yt  is  a function  of  f only*  It  may  be  denoted  by  H(  £ ).  This 

function  fulfills  the  differential  equation 


jrfli 


(23) 


and  may  be  defined  by  the  boundary  conditions 

H « 1 for  £ -*>  - oo 

H = 0 for  ^ » 1 (23a) 

The  function  H(  ^ ) can  be  expressed  by  a linear  combination  of  the 
symmetric  and  antisymmetric  solution  G for  A ■ 1/16,  The  symmetric 
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solution  is  simply  G ■ 1.  Since  Ga(J  , l/l6)  is  zero  for 
H is  expressed  by  J 


H '/+  's-G'a,  A) 


co, 


here  is  a coefficient  which  is  determined  by  the  condition  H * 0 
at  ^ « 1*  Using  Equation  (13)  one  obtains  for  the  value  of  H at  £ ■ 1 
according  to  the  last  equation 


H-i  * 


ci 


rfi) 


r(ii-rfi) 


Hence 


5 


z dtl 

/?  ■/?  di) 


and 


H **  / ~ 


£ ciii  /*  //  -i  j 


m 


To  the  expression  H a superposition  of  the  solutions 

and  Yu.  with  undetermined  coefficients  will  be  added*  Parti- 
cular solutions  will  not  occur  since  they  would  introduce  a 

singularity  at  point  0*  The  particular  solutions  and  A* 

-fa 

are  also  singular  at  point  0,  because  of  the  factor  g which 

occurs  (see  Equation  (9))j  however,  the  final  solution  will  not  show 
this  singularity*  Thus  the  representation  of  the  inner  solution  will 
be  given  by 
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¥*£%(£)  • & ft,  f'i+'it’) 

j „ 

/ *° 

+(fj  JftW'G ft, 

(25) 

The  Constants  su  and  the  functions  oC  ( ? ) and  /3  ( >*  ) will  be  deter- 
mined by  the  matching  process  along  ■ £e. 

For  the  representation  of  the  outer  solutions  the  boundary  condi- 
tions at  the  wedge  surface  BCD  must  be  taken  into  account.  One  first 
constructs  a family  of  particular  solutions  which  fulfills  along  the 
line  OABCD  the  boundary  condition  * 0 and  which  have  a singularity 
at  point  0 given  by  U',  . These  particular  solutions  have  the  form 


(26) 

where  the  functions  fulfill  the  conditions  of  Tricani's  uniqueness 

and  existence  proof  (Reference  3)*  It  has  been  shovm  in  Reference  2 
that  any  solution,  which  fulfills  outside  of  a line  = ? the 

boundary  condition  ■ 0 along  ABCD,  can  be  represented  by  a super- 
position of  the  expressions  (26).  Thus  the  outer  solution  is 
represented  by 

?-£**'*  '% 

(27) 

The  coefficients  b^  will  be  determined  by  the  matching  process.  If 
the  choking  Mach  number  is  quite  close  to  one,  the  line  ° q Q 

lies  rather  closely  to  the  origin  and  the  singular  terms  in  the 
expressions  (26)  will  prevail,  so  that  for  the  purpose  of  matching 
along  the  line  o m P Q the  outer  solution  can  be  replaced  by  the 
expression 
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(27a) 


00  -/V  +*! 


By  matching  the  inner  and  outer  representations  for  the  function  y. > 
and  their  derivatives  with  respect  to  f along  the  line  e>  = f 
one  obtains  the  equations 


GO 

H({  ] +L  % G ff,  A i *3>>fj  +J°c(>j  (3  7/,  - 'fi'VoV 

>*-/ 


(2&l) 


oo  » 

27  <?  7f,  f-i  *m/j  -sty {-,{«»>-* 

" *27  /■/ ^ 7)7  TV  +1/>Yy 

* * 


(22b) 

The  function  H can  be  developed  in  terms  of  the  functions  G by  means 
of  Equation  (22).  By  suitable  integration  by  parts  which  utilize  the 
differential  equation  (23)  for  H and  the  differential  equation  (7)  for 
the  G and  furthermore  the  boundary  conditions  for  G and  H (Equation 
(23a))  one  finds 


Oo  " 
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and 


/ 


With  Equations  (lit.)  and  (16)  one  thus  finds  as  the  representation  for 
H 


H 


_jf  £j 'JA.> 


A-/ 


/V4K 


J .Jr 

J ,{)  At?) 

o 


Inserting  this  expression  into  Equation  (27a)  one  finds  by  equating 
the  coefficients  of  the  different  expressions  G 


**  ‘ *4 -</•<*  ' * 

c£  . _£  — ' 

Similarly  from  Equation  (27b) 

a,  (-$  +h)  - ~ +hJ-  &/, 

n 

~d  • ° 


Hence 


ah  . — 


/-i  *A)_ 


^ ^ ~~tU‘  ^ {'£ 


' 


oC(pJ  J 


^ ^/4?/  x P 


AV  - ; — • * o/ 
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This  determines  the  solution.  We  are  particularly  interested  in  the 
representation  of  the  inner  region,  since  it  gives  the  pressure  dis- 
tribution along  the  w±nd-4>unnel  wall.  It  is  given  by 


SECTION  VI 

TRANSFORMATION  OF  EQUATION  (29) 


A direct  evaluation  of  the  last  expression  would  be  quite 
cumbersome  because  of  the  occurrence  of  the  complex  function  A(  )>  ), 

defined  in  Equation  (IS),  One  observes,  however,  that  the  integrand 
is  an  analytic  function  of  9 and  therefore  by  a deformation  of  the 
path  of  integration  in  the  complex  V -plane  one  can  eventually  obtain 
a more  convenient  expression.  Actually,  the  solution  will  finally 
appear  as  an  infinite  series  which  also  displays  more  clearly  the 
behavior  of  the  solution  near  point  0,  This  transformation  will  be 
carried  out  in  the  following  solutions  but  instead  of  the  expression  V/ 
we  shall  determine  directly  the  x coordinate  along  the  wind-tunnel  wall 
as  a function  of  £ , for  this  is  the  quantity  in  which  we  are  mainly 
interested. 

In  the  present  simplification  of  the  hodograph  equation  one  has 
the  relation 

- (/+•>*  f ft 
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for  simplicity  the  stream  density  at  sonic  "speed  (i.e*,  the  flux  per 
unit  area)  has  been  put  equal  to  1*  The  functions  behave  along  the 
line  ^ ■ - oo  as  (-  f )**3/2. 

With 


If  - A 

and  If  - 0 

20 


-J/. 

for  / 


00 


one  thus  finds  from  Equation  (29).  inserting  at  the  same  time  the 
expression  for  H from  Equation  (24), 


i 


- z rty 

/S'  /T  rtf! 


O© 

T±  —L 

£-•2  /i 


(-£ 


- f m 


h*f 


A ~iJ<  H+M) 


+ / 1 SJ*  > /%  + ar  £/ 


Here  £ is  the  value  of  \ for  the  point  A,  ^ is  negative. 
Then  x for  the  wind-tunnel  wall  is  found  by  an  integration;  x = 0 
may  lie  at  the  intersection  of  the  sonic  line  with  the  wall. 


X - 


_ / 


A ~/t ) dt  f* 
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The  integral  -which  occurs  on  the  right  hand  side  may  now  be  trans- 
formed* It  may  be  denoted  by  J. 


_ f j. 


(31) 


Inserting  into  the  last  expression  the  value  of  A from  Equation  (IS), 
one  finds 


J. 


**J  -p  rfu?)  • r(-zc») 

d 


Because  of  the  relations 


jcn'fZ'fi'V  + 7j-)h 1,9  J * / co-sttoP 

rfr>  +£)r/-<:>+£)  , /yt- , . -7  ±j 

V*-  + >//*?  ¥ 


r(zcyj>  rf-ioy)  - 


this  expression  cai  be  simplified  to 


J*  -jp.-tJ' rf~iv  ~ii ) rfcV  +£) • 


One  can  write  it  as 


± ReJ^  (z^v) ' rf*  i ft*  +*.)'  + i)  l Jy 


IS 
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We  set 


i V « u,  then 


too  2/ 

/ YXu-  r(u'i) 


The  path  of  integration  is  along  the  positive  imaginary  u axis. 
Further  simplifications  are  carried  out  by  means  of  the  formula 


ru>  rf,-z) , 


K 

SvnXZ. 


and  the  foil  wring  relations  for  trigonometric  functions 


ClH  Zl4_-  Z StAv( U + . St iv  f-K,  + 

JWoC  -Ah/«C  +g)  . SV*v( aC  <-f)  , ^ 


Thus  one  obtains  finally 


c °o 

£ Ct 


V(  y ic  y 


ChzALOk&L  AJk_  ttj 


We  shall  deform  the  path  of  integration  from  the  imaginary  axis  to  the 
positive  real  axis.  Obviously  all  poles  of  the  integrand  lie  at  the 
real  axis.  This  deformation  will  be  possible  if  the  contribution  to 
the  integral  along  a line  which  leads  at  a far  distance  from  the 
imaginary  axis  in  the  complex  u-plane  to  the  real  axis,  drops  out. 

By  means  of  the  asymptotic  expression  for  the  /"*  -function 
which  are  valid  in  the  entire  u-plane  with  the  exception  of  the  nega- 
tive real  axis,  one  finds  that  the  combined  P -functions  in  the 
bracket  of  the  last  expression  give  a behavior  as  u“l*  The  combined 
trigonometric  functions  have  a periodicity  2 with  respect  to  the  real 
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part  of  u,  -with  respect  to  the  imaginary  part  of  u it  behaves  as 
e-  2 \ Since  £</  the  expression  S'H.  decreases  expon- 


entially  in  the  positive  real  u-direction.  On  the  basis  of  these 
properties  it  is  not  difficult  to  find  a suitable  family  of  contours 
connecting  the  real  to  the  imaginary  axis  and  avoiding  the  poles  which 
lie  on  the  real  axis,  such  that  the  contribution  of  the  integral  tends 
to  zero  if  u along  the  contours  exceeds  any  limit.  Thus  the  integra- 
tion can  be  carried  along  the  path  in  the  u-plane  shewn  in  Figure  5# 
where  the  singular  points  are  excluded  by  small  circles.  Now  because 
of  the  factor  l/i  in  front  of  the  integral  the  portions  of  the  path 
of  integration  along  the  real  axis  do  not  contribute  to  the  real  part 
of  the  last  expression  and  only  the  contribution  of  the  circles  around 
the  poles  is  of  importance.  Thus  one  finds 


J * • fouduts  of- 

KJ  ' 


r(u->g)r(H+& 


(32) 


The  poles  of  the  integrand  along  the  positive  real  axis  lie  at 

u = 1/12,  u«-l/U+h(h-l,  2...  ) 
u - - 1/12  + h/3  (h  - 1,  2 ...  ) 

For  the  points  u = -1/12  + h/3  one  finds  by  direct  evaluation 


/ ) 

s*»,/u+%) 


Thus  the  residues  at  the  point  u = -1/12  + h/3  are  found  to  be 


-r(j-i)  r[§H)  /A*-* 

jrrfp  o rff'J)  V 


/}  ° / 2 < >< 
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For  the  residues  at  points  u *»  - l/b  + *h  one  obtains 


r(f-fh)  r(i-fh/  hJ 

Finally  the  residues  at  the  point  u = 1/12  is  found  to  be 


{W 

ff-na-rui 

Inserting  into  Equation  (30),  the  Equations  (31) , (32)  and  the 
expressions  for  the  residues  just  found,  introducing  furthermore  for 
the  the  expressions  (14),  and  for  the  the  expressions  (15), 

one  finds  that  the  residue  at  the  point  u «*  l/l2  cancels  the  contri- 
bution to  x caused  by  the  function  H (this  is  the  first  term  of  the 
bracket  in  Equation  (30)  and  that  the  residues  at  the  points  u ® 

- 1/lj.  + h cancels  the  contribution  of  the  sum  in  the  bracket  of 
Equation  (30).  Thus  the  result  assumes  the  form 


jr 


oo 


ntiL  1 rA  *JL 

-/YM) 


h 


(33) 


SECTION  VII 
DISCUSSION 


The  expression  for  has  been  found  with  the  assumption  that 
^ «*  1 along  the  wind-tunnel  wall.  For  a wind  tunnel  of  height  H, 
the  result  must  be  multiplied  with  H/2.  is  the  value  of  £ at 

point  A and  point  A corresponds  to  the  choking  Mach  number.  By  the 
definition  of  ^ , Equation  (1)  and  the  relation 
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one  finds 


M-t 


Vj 

i (/->')  ■? 


Thus  Equation  (33)  can  be  written  as 


H it 


/-  /7 


) 


where 


/ rd-y  rtf+i)  //-//  f 
£ rtf*,). rtf  tf)  I t-Huu/ 


The  function  E is  shovm  in  Figure  6.  One  will  surmise  that  the  con- 
tinuation of  the  pressure  distribution  into  the  supersonic  region  is 
the  analytical  continuation  of  the  above  expression.  To  show  this, 
one  must  follow  the  solution  through  the  supersonic  legion  of  the 
hodograph.  One  will  remember  that  the  flow  in  the  vicinity  of  the 
centerline  of  a Laval  nozzle  maps  into  three  sheets  of  the  hodograph, 
and  that  the  flow  at  the  wind-tunnel  wall  is  equivalent  to  such  a flow. 


In  the  transonic  approximation  it  is  possible  to  interpret  the 


expression 


1 - M 


1 - 


as 


P - P* 


choking 


^choking  " * 


where  p is  the  pressure 


ptt  is  the  pressure  at  sonic  velocity.  The  investigations  have  been 
carried  out  for  a wedge  in  a wind  tunnel,  but  they  have  a more  general 
validity*  The  shape  of  the  body  is  expressed  by  the  functions 
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which  occur  in  Equation  (26).  However,  by  the  assumption  of  a very 
wide  tunnel  this  function  is  eliminated  from  the  future  computations. 
So  the  result  obtained  in  Equation  (33)  is  valid  for  any  body  shape, 
if  the  wind  tunnel  is  sufficiently  wide. 

In  this  representation  the  governing  parameters  are  the  wind- 
tunnel  width  and  the  choking  Mach  number.  In  the  paper  of  Marschner, 
Reference  1,  a connection  has  been  established  between  the  choking 
Mach  number  and  the  dimensions  of  a diamond  shaped  airfoil 


/- 


//  I*/*-  *//F 

,uUhI  4 


Thus  the  wall-pressure  distribution  can  be  connected  to  the  size  of  a 
diamond  profile. 

In  Marschner fs  paper  a criterion  for  the  validity  of  the  approxi- 
mation has  been  derived.  It  is  given  by  the  fact  that  the  value  of  p 
must  be  small  in  comparison  to  the  value  of  f at  the  point  C,  i.e., 
at  the  sonic  point  of  the  shoulder  characteristic.  This  criterion 
requires  that 


L/s  - Vr 

HI  0-  «' 


Thus  the  accuracy  of  the  results  depends  not  only  upon  the  choking 
Mach  number  but  also  upon  the  slenderness  of  the  profile.  Also  from 
a mathematical  point  of  view  the  result  offers  some  interest.  In 
contrast  to  Equation  (29),  only  integral  powers  of  £/»£<>  occur  in 

Equation  (33)*  Thus  the  singularities  at  point  0 introduced  in  the 
individual  particular  solutions  Ykh  and  fwh  cancel  out  by  the 

superposition  in  the  Fourier  integral. 

If  one  has  to  investigate  the  vicinity  of  point  0 one  will  ask 
which  particular  solutions  of  the  form  (5)  would  be  suitable  to 
represent  the  solution  in  the  vicinity  of  point  0 if  one  considers 
the  boundary  condition  imposed  by  the  nature  of  the  problem.  (In  this 
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regard  the  conditions  imposed  in  the  formulation  of  the  eigen  value 
problem  are  quite  artificial.)  These  conditions  require  that  the 
function  is  antisymmetric  ■with  respect  to  AO  (except  for  a solution 
\t>  » 1)  and  furthermore  we  know  that  no  singularity  propagates  along 

tne  limiting  characteristic  DO.  Finally  point  0 must  not  map  to 
infinity.  By  the  condition  « 0 along  OA  we  are  restricted,  as 
previously,  to  solutions  . Their  behavior  along  the  limiting 
characteristic  DO  can  be  investigated  by  writing  the  arising  expressions 
y in  the  following  form  which  can  be  obtained  by  a combination  of 
Equations  (5)  and  (13). 


ff 


■ G ({ A) -i  ■ rit)lrf$+£)  r($+i I 


The  hype rgeome trie  functions  which  occur  in  Equation  (13)  have  been 
replaced  by  1 since  ^ is  close  to  1.  The  last  equation  shows  that 

the  limiting  Mach  wave  would  map  to  infinity  (for  ^ < 0),  or  a 

singularity  would  propagate  along  it  (for  /X  > 0)  if  the  second  term 
in  the  bracket  is  different  from  zero.  From  the  condition  that  the 
coefficient  of  the  second  term  vanishes  one  obtains  certain  values  of 
y (incidentally,  they  are  different  from  the  values  of  A for  the 
eigen  functions).  If  one  determines  for  the  particular  solutions 
determined  with  these  values  of  , the  expressions  for  x along  the 
line  OA  one  finds  the  powers  of  encountered  in  Equation  (33 )«  So 

the  form  of  Equation  (33)  is  not  unexpected.  Nevertheless,  it  does  not 
seem  practical  to  introduce  the  particular  solution  which  appears  in 
Equation  (33)  right  at  the  beginning  for  then  it  would  have  been  very 
difficult  to  carry  out  the  matching  process  along  the  line  p ■'  <?  Q. 

From  the  form  (33)  one  can  naturally  obtain  a representation  for 
y along  the  limiting  Mach  wave  DO.  It  proceeds  also  in  power  of  *l0 


and  shows  that  is  different  from  zero  along  the  limiting  Mach  wave. 
Thus  the  artificial  boundary  condition  y,  = 0 along  £ - Cg  will  not 

cause  the  solution  Y'  to  be  zero  along  DO  after  the  limiting  process 
c,,  -*■  1 has  been  carried  out. 
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2t  Hodogrsph  Mop  of  Figuro  1 


Figaro  3*  Hop  of  Figaro  1 in  tho  7 , 6 -flono 
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WIND-TUNNEL  WALL 


Figure  J*r  Body  With  Rear  Part  Which  Does  Not  Reflect  Waves 
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IMAGINARY  AXIS 
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Vigor*  5t  Path  of  Integral  In  Cmplax  u-Plana 


Figure  6t  Mach  Number  Distribution  Along  Tind  Tunnel  lfcll 
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